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Abstract 

Inclusive b(1020)-meson production in neutral current deep inelastic e'^p scat¬ 
tering has been measured with the ZEUS detector at HERA using an in¬ 
tegrated luminosity of 45 pb“^. The 0 mesons were studied in the range 
10 < < 100 GeV^, where is the virtuality of the exchanged photon, and 

in restricted kinematic regions in the transverse momentum, px, pseudorapidity, 
p, and the scaled momentum in the Breit frame, Xp. Monte Carlo models with 
the strangeness-suppression factor as determined by analyses of e'''e“ annihila¬ 
tion events overestimate the cross sections. A smaller value of the strangeness- 
suppression factor reduces the predicted cross sections, but fails to reproduce the 
shapes of the measured differential cross sections. High-momentum 0 mesons in 
the current region of the Breit frame give the first direct evidence for the strange 
sea in the proton at low x. 
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1 Introduction 


The total quark content of the proton has been well determined iai§ through analyses 
of inclusive deep inelastic scattering (DIS) data. However, the flavour decomposition of 
the sea is less well known. So far, experimental constraints on the strange-quark content 
of the nucleon have come from fixed-target neutrino experiments 0 ] , which indicate that 
the ss is suppressed with respect to the uu and dd sea by a factor of about two. This paper 
reports a study of the production of 0-mesons in neutral current e'^p DIS and explores 
its sensitivity to the strange sea of the proton at low x. 


Several mechanisms lead to 0-meson production in DIS. The 0 meson, which is a nearly 
pure ss state, can be produced by the hadronisation of a strange quark created in the 
hard scattering process of a virtual photon on the strange sea of the proton, 7*5 —>• 
s, as illustrated in Fig. |]a). The underlying hard-scattering process is either zeroth 
order in QCD, namely the quark-parton model (QPM), or hrst order, 7*5 —> sg, the 
QCD Compton reaction (QCDC). Another source of strange quarks is boson-gluon fusion 
(BGF), '-y*g —> ss. Fig. |J)). In contrast to the QPM and QCDC processes, the rate of 
BGF events is related to the density of gluons in the proton and is, therefore, not directly 
dependent on the intrinsic sea-quark content of the proton. The hadronisation process 
alone, without strange quarks being involved in the hard scattering, contributes to the 
production of 0 mesons, as shown in Figs. |^c)-d). In this case, 0 mesons are formed from 
strange quarks created during hadronisation. Hadronisation of strange quarks produced 
in higher-order QCD reactions related to the splitting of gluons, g ss, and the decay 
of higher-mass states, such as the Dg meson (Fig. He)), also contribute. In addition, 
diffractive scattering can produce 0 mesons in the final state (Fig. |l|f)). 


Strange-particle production in inclusive DIS has been studied at HERA using mesons 
and A baryons BO. However, their production rates are dominated by the fragmenta¬ 
tion process and by the decays of high-mass states, and are, therefore, insensitive to the 
presence of strange quarks in the hard scattering process. For 0 mesons, the contribution 
from resonance decays is relatively small. Furthermore, selecting 0 mesons with large 
longitudinal momenta in the Breit frame |T^ enhances the contribution from the QPM 
process of Fig. |T]a). 


In this study, the 0 mesons were identihed through the decay 0 — K^K~. Their differ¬ 
ential cross sections are presented as functions of = —q^ = —{k — k'Y and Bjorken 
X = Q‘^/{2Pq), where k and k' are the four-momenta of the initial and scattered lep¬ 
ton and P is the four-momentum of the incoming proton, as well as other variables that 
characterise the 0 -meson production. 
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2 Properties of 0 mesons in the Breit frame 


The Breit frame [^] provides a natural system to separate the radiation of the outgoing 
struck quark from the proton remnant. In this frame, the exchanged virtual boson with 
virtuality Q is space-like and has a momentum q = {qo, qxsy Qzb) = (0,0,0,—Q). 
In the QPM, the incident quark has pzj^ = Q/‘2. and the outgoing struck quark carries 
Pzb — All particles with negative pz^ form the current region. These particles are 

produced by the fragmentation of the struck quark, so that this region is analogous to a 
single hemisphere of an annihilation event. 

The (^-meson cross sections are presented as a function of the scaled momentum, Xp = 
2p/Q, where p is the absolute momentum of the 0 meson in the Breit frame. In the QPM 
process, 7*5 —>■ s, this variable is equal to unity for the s-quarks in the current region. As 
a consequence, leading 0 mesons in the current region with Xp values close to unity are 
a measure of the hard scattering of a virtual photon on the strange sea. Gluon radiation 
and the fragmentation process generally lead to particles with Xp < 1 , and, much less 
frequently, to Xp > 1 . 

In the target region, Xp can be signihcantly larger than unity. This is because the 
maximum momentum of the proton remnant in the QPM is Q{1 — x)/2x, therefore 
^max _ (^i — x)/x. The 0 mesous in the target region are mostly produced by the hadro- 
nisation processes of Fig.lc)-d), as well as the hadronisation of strange quarks from the 
BGF diagram of Fig. lb). 


3 Data sample and analysis procedure 


3.1 Experimental setup 


During the 1995-1997 period, 45.0 ± 0.7 pb ^ of data were taken with the ZEUS detector 
with a positron beam energy of 27.5 GeV and a proton beam energy of 820 GeV. 


ZEUS is a multipurpose detector described in detail elsewhere |^. Of particular impor¬ 
tance in the present study are the central tracking detector and the calorimeter. 


The central tracking detector (CTD) is a cylindrical drift chamber with nine super¬ 
layers covering the polar-angle|^ region 15° < 9 < 164° and the radial range 18.2 — 79.4 cm. 


^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the “forward direction”, and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — ln(tan |), where the polar angle, 0, is measured with respect to the proton beam 
direction. 
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Each superlayer consists of eight sense-wire layers. The transverse-momentum resolution 
for charged tracks traversing all CTD layers is a{pT)/pT = O.OOSSpT©0.0065©O.OOld/pr, 
with pt in GeV. 


The CTD is surrounded by the uranium-scintillator calorimeter, CAL |^, which is di¬ 
vided into three parts: forward, barrel and rear. The calorimeter is longitudinally seg¬ 
mented into electromagnetic and hadronic sections. The smallest subdivision of the CAL 
is called a cell. The energy resolution of the calorimeter under test-beam conditions is 
(JeIE = 0.18/a/E for electrons and (TeIE = 0.35/a/E for hadrons (with E in GeV). 


The position of positrons scattered at small angles to the positron beam direction was 
measured using the small-angle rear tracking detector (SRTD) |^,^]. The energy of 
the scattered positrons was corrected for the energy loss in the material between the 
interaction point and the calorimeter using a presampler (PRES) 


3.2 Kinematic reconstruction and event selection 


The scattered-positron candidate was identihed from the pattern of energy deposits in 
The kinematic variables, and x, were reconstructed by the following 


the CAL 
methods: 
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the electron method (this method is denoted by the subscript e) uses measurements 
of the energy and angle of the scattered positron; 


the double angle (DA) method relies on the angles of the scattered positron and 
of the hadronic energy flow; 


• the Jacquet-Blondel (JB) method is based entirely on measurements of the hadronic 
system. 

The DIS event selection was based on the following requirements: 


• Ep > 10 GeV, where Ep is the energy of the scattered positron in the calorimeter 
after the correction by the PRES; 

• 10 < Qg < 100 GeV^. The upper cut on Ql was used to reduce the combinatorial 
background in the 0-meson reconstruction; 

• 40 < 5 < 60 GeV, where 6 = Ei{l — cos9i), Ei is the energy of the ith calorimeter 
cell, 9i is its angle, and the sum runs over all cells. This cut further reduces the 
background from photoproduction and events with large initial-state radiation; 

•2/e A 0.95, to remove events with fake scattered positrons; 

• 2/jb > 0.04, to improve the accuracy of the DA reconstruction used in systematic 
checks; 
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• a primary vertex position, determined from the tracks fitted to the vertex, in the range 
I ^vertex |< 50 cm, to reduce backgronnd events from non-ep interactions; 

• the impact point {X, Y) of the scattered positron in the calorimeter must be within 
a radius y/X^~+Y^ >25 cm. 

The reconstruction of the Breit frame and the and x variables was performed using 
the electron method, since it has the best resolution at the relatively low values of this 
data set. 


4 Selection of <j) candidates 


Charged tracks measured by the CTD and assigned to the primary event vertex were 
selected. Tracks were required to pass through at least three CTD superlayers and have 
transverse momenta pr > 200 MeV in the laboratory frame, thus restricting the study to 
a CTD region where track acceptance and resolution are high. 

All pairs of oppositely charged tracks were combined to form the 0 candidates. The 
tracks were assigned the mass of a charged kaon when calculating the invariant mass, 
M{K~^K~), of each track pair. The events with 0-meson candidates were selected using 
the following requirements: 

• 0.99 < M{K+K-) < 1.06 GeV; 

• > 1.7 GeV and —1.7 < < 1.6, where p^ is the transverse momentum and is 

the pseudorapidity of the 0 meson in the laboratory frame. 

The asymmetric cut on was used to avoid the very forward region that has large track 
multiplicities, resulting in high combinatorial backgrounds. 

Figure shows the invariant-mass distribution for 0 candidates in the range 10 < < 

100 GeV^. The invariant mass for the leading 0 mesons in the current region of the 
Breit frame, 0.8 < Xp < 1.1, is presented in Fig. |b). For the latter case, the DIS events 
containing the 0-meson candidates have x < 0.006. The solid line in each hgure is a £t 
using a relativistic Breit-Wigner (BW) function convoluted with a Gaussian distribution 
plus a term describing the background: 

F{M) = {BW) ® (Gaussian) + a{M — 2mx)^, 


where a and b are free parameters and is the kaon mass. The fit function contains hve 
free parameters: normalisation, peak position, width of the Gaussian distribution, and two 
parameters describing the background. When the peak position was left free, the resulting 
£t gave 1019.2±0.3 MeV, in agreement with the PDG value of 1019.456±0.020 MeV p9| . 
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The width of the Gaussian was 1.6 ± 0.3 MeV, consistent with the tracking resolution. In 
order to improve the stability of the £t for the calculations of the differential cross sections, 


the mass peak and width of the Breit-Wigner function were hxed at the PDG values 


The total number of 0-meson candidates determined from this £t was 4950 ± 214, while 
the number of 0-meson candidates for the high Xp region was 181 ± 28. 


5 Event simulations 


A good understanding of hadronisation is a pre-requisite for the interpretation of the mea¬ 
sured inclusive 0-meson cross sections. At present, only Monte Garlo (MG) models based 
on leading-order QGD are available to compare with the experimental results, so that the 
predictions for the rates of ss production are plagued by large model-dependent uncer¬ 
tainties. In MG models based on the Lund string fragmentation , the production ratio 
of strange to light non-strange quarks is parameterised by the strangeness-suppression 
factor. As = Ps/Pu,d, where Pg {Pu,d) is the probability of creating s {u,d) quarks in the 
colour held during fragmentation. The processes shown in Figs. |^) and b) are propor¬ 
tional to As, while the contributions illustrated in Figs, ^c) and d) are proportional to 

In annihilation, the production of 0-mesons has been well described using As = 0.3 
3T||. However, there are new indications that a larger value. As — 0.4, may be needed ||35|| , 


or even that a single value cannot accommodate all of the SLD strangeness-production 
data [^]. When using the same hadronisation model in e'^p scattering, the measured 
and A production rates in DIS [|T^|T^ and photoproduction indicate the need for a 
smaller value. As — 0.2. 

The measured cross sections were compared to various leading-order MG models based 
on the QGD parton-cascade approach, to incorporate higher-order QGD effects, followed 
by fragmentation into hadrons. The MG events were generated with LEPTO 6.5 ||38|| , 
ARIADNE 4.07 and HERWIG 6.2 using the default parameters in each case. The 
fragmentation in LEPTO and ARIADNE is simulated using the Lund string model [OT as 


implemented in PYTHIA [41|, whereas the hadronisation stage in HERWIG is described 


by a cluster fragmentation model 


The acceptance was calculated using ARIADNE, which was interfaced with HERAGLES 
4.5.2 |0| using the DJANGOH program |45] in order to incorporate first-order electroweak 
corrections. The generated events were then passed through a full simulation of the 
detector using GEANT 3.13 [Q and processed with the same reconstruction program as 
used for the data. The detector-level MG samples were then selected in the same way as 
the data. 
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The natural width of the Breit-Wigner distribution for (/)-meson decays was set to the 
default value in LEPTO and ARIADNE. The HERWIG model sets the particle- 
decay width to zero and is therefore less realistic for the acceptance calculations. The 
HERWIG model was used only for comparisons with the hnal cross sections. 

The inclusive 0-meson sample contains a contribution from diffractive processes, which is 
not well simulated in the MG models mentioned above. These processes are characterised 
by a rapidity gap, chosen as < 2, where is dehned as the pseudorapidity of the 
energy deposit in the GAL above 400 MeV closest to the proton direction, and by the 
presence in the GTD of only a few tracks. Diffractive events with 0 mesons were generated 
with PYTHIA 5.7 and passed through the same simulation of the detector as for 
inclusive MG events. The MG distributions were £t to the data by varying the fraction 
of the diffractive 0-meson events from PYTHIA and minimising the to obtain good 
agreement for the multiplicity of charged tracks in the GTD. The fraction of PYTHIA 
events needed to obtain good agreement between data and MG was 2.7 ± 0.2% of the 
total number of reconstructed 0-meson events. It was verihed that this fraction gives a 
satisfactory description of the 0- meson events for r^max < 2. 


6 Definition of cross sections and systematic uncer¬ 
tainties 


The 0-meson cross sections were measured in the kinematic region 10 < < 100 GeV^, 

2 ■ 10“^ < X < 10“^, 1.7 < < 7 GeV and —1.7 < < 1.6. The cross sections as a 

function of a given observable, Y, were determined from 


da N 

1y ^ A-C-B-AY ’ 


where N is the number of events with a 0-meson candidate in a bin of size AY, A is the 
correction factor (which takes into account migrations, efficiencies and radiative effects 
for that bin) and C is the integrated luminosity. The branching ratio, B, for the decay 
channel 0 —> K^K~ was taken to be 0.49210™® 


The acceptance for each kinematic bin was calculated as where (.4,®®“) is 

the reconstructed (generated) number of events with 0 mesons. For the calculation of the 
acceptance, 2.7% of the total number of inclusive DIS events generated with ARIADNE 
were replaced by diffractive events from PYTHIA. While the contribution from diffractive 
0 -meson events is negligible for the full phase-space region, it is important for the high 
Xp region in the Breit frame, since 72% of the diffractive 0-meson events have Xp > 0.8. 
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The systematic uncertainties of the measured cross sections were estimated from the 
following (the typical contribution from each item to the uncertainty of the total cross 
section is indicated in parentheses): 

• event reconstruction and selection. Systematic checks were performed by changing the 
cuts on He, yjB, ^ and the vertex-position requirement: ye < 0.90 (—0.1%), yjB > 0.05 
(—0.05%), 42 < 5 < 58 GeV (—0.3%), | ^vertex |< 45 cm (-1-0.4%). The radius 
cut for the position of the scattered positron in the calorimeter was raised by 1 cm 
(—0.5%). The minimum accepted energy of the scattered positron was increased by 
1 GeV (—0.1%). The positron energy scale was changed within its ±2% uncertainty 

ao:7%); 

• the DA method was used to reconstruct the Breit frame (-1-0.3%) and the kinematic 
variables (-1-0.08%); 

• the minimum transverse momentum for A'-meson candidates was raised by 100 MeV 
(-1-0.6%). Tracks were required to have \ rj \< 1.75, in addition to the requirement of 
three GTD superlayers (-1-0.02%); 

• the form of the background in the hts was changed to a second-order polynomial 
function (-|-0.4%); 

• the fraction of diffractive ())-meson events in the Monte Garlo sample was varied in the 
range 1.9 - 3.5% (±0.03%). 

The overall systematic uncertainty for the differential cross sections was determined by 
adding the above uncertainties in quadrature. The normalisation uncertainty due to that 
of the luminosity measurement, which is 1.6%, was only added to the overall systematic 
uncertainty for the total 0-meson cross section. The uncertainty in the 0 —>■ K^K~ decay 
branching ratio was not included. 


7 Results 

The overall 0-meson acceptance for 10 < < 100 GeV^, 2 ■ 10“^ < x < 10“^, 1.7 < < 

7 GeV and —1.7 < estimated with DJANGOH, was 45%. The total 0-meson 

cross section in this region is 

a{e^p —>• e^(f)X) = 0.507 ± 0.022(stat.)lQ;QQg(syst.) nb. 

This cross section is lower than that predicted by the LEPTO (0.680 nb) and ARIADNE 
(0.701 nb) models with the GTEQ5D structure function and with the LEP default value of 
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the strangeness-suppression factor, = 0.3. The HERWIG 6.2 model for neutral current 
DIS processes underestimates the measured 0-meson cross section, predicting 0.36 nb. 

In previous studies of neutral kaons and A baryons at HERA 0,00, it was found 
that decreasing from its standard value improved the agreement between the Lund MC 
models and the data. A smaller value of the strangeness-suppression factor, A^ = 0.22, 
resulted in an inclusive 0-meson cross section of 0.501(0.509) nb for LEPTO (ARIADNE), 
which agrees well with the present measurement. Therefore, A^ = 0.22 was used as the 
default for LEPTO and ARIADNE in the following comparisons. A comparison of the 
data with the predicted cross sections gave an uncertainty of ±0.02 on the A^ value used 
in this analysis. 


7.1 Differential 0-meson cross sections 


Figure ^ shows the differential cross sections as a function of Xp and compared 

to the LEPTO, ARIADNE and HERWIG models using the GTEQ5D parton distribution 
functions^. The measured cross sections are compiled in Table |I|a) - e). The Xp cross 
sections are shown separately for the current and the target regions. The 0-meson cross 
sections in the current and the target regions of the Breit frame are distinctly different: 
the data are concentrated at Xp around ~ 0.5 in the current region, and at ~ 1 in the 
target region. 


The MG models based on the Lund fragmentation with A^ = 0.22 reasonably well repro¬ 
duce the pi^ and distributions. Signihcant differences exist for the distributions of 
in the laboratory frame and Xp in the current region of the Breit frame. In the target 
region, the MG models underestimate the cross sections. If A^ = 0.3 is used, the MG 
models based on the string fragmentation agree well with the data in the target region, 
but signihcantly overestimate the cross sections in the current region. 


In addition to varying the A^ values, different methods to tune the Lund MG models were 
considered, all of which had a negligible effect on the LEPTO predictions. In particular, 
the contribution to the 0 cross section from charm events, mainly due to D and decays, 
was investigated using AROMA This model produces charm quarks exclusively 

through the BGF mechanism, and reproduces the measured D*^ cross sections in DIS 


According to AROMA, charm decays account for 20% of the 0 mesons, contributing 
mainly in the target hemisphere. This fraction is larger than that predicted by LEPTO, 
but it is not sufficient to explain the observed discrepancies. For leading 0 mesons {xp > 
0 .8) in the current region, charm events give a negligible contribution. 


^ The leading-order set, CTEQ5L, gave very similar results. 





In order to disentangle the different contributions to the 0-meson production and to in¬ 
vestigate the observed discrepancies, the MC samples were divided into a few subsamples. 
Figure illustrates the contributions of QPM/QCDC interactions on an s or s quark of 
the proton sea. In this case, a struck s or s quark produces a 0 meson after the hadronisa- 
tion process. The 0-meson cross section in the current region of the Breit frame contains 
a significant fraction of events produced by hard scatterings of the virtual photon on the 
strange sea. This fraction rises with increasing and Xp values, while the contribution 
to 0-meson production from strange quarks produced solely in the hadronisation process 
becomes negligible for Xp > 0.8. In contrast, the target region contains a small contribu¬ 
tion from the QPM/QCDC events, since the second s or s from an ss pair participating 
in the interaction usually escapes undetected in the very forward region. 

Figure also indicates the contribution of BGF processes in which the flavour of the 
produced quark is s or s. The fraction of these BGF events is larger in the target region 
than in the current region. 

7.2 The 0)-meson cross section as a function of x 

Production of 0 mesons was investigated as a function of x. The s-quark density increases 
with decreasing x; however, the BGF contribution also increases with decreasing x due to 
the rise of the gluon density. Thus, the 0-meson cross section as a function of x depends 
on both the strange sea and the gluon density. 

The differential cross sections as a function of x for two regions, 10 < < 35 GeV^ 

and 35 < < 100 GeV^, are shown in Fig. Table gives the values of the cross 

sections. The 0-meson differential cross section increases as x decreases down to the 
kinematic limit. The LEPTO and HERWIG MGs reproduce this rise. The LEPTO model 
shows the contributions of events in which a 0 meson is produced after hadronisation of an 
s (s) quark emerging from the hard interaction. The contributions from the QPM/QGDG 
and BGF processes rise with decreasing x due to the rise of the s-quark and the gluon 
density in the proton. 

7.3 Leading 0 mesons 

The MG predictions for leading 0 mesons {xp > 0.8), usually corresponding to high in 
the laboratory frame, have small uncertainties both in the simulation of the QGD processes 
and in the hadronisation mechanism; for a given strange-sea density, the scattering of the 
virtual photon on a strange quark is described by the QED process, 7 *s —> s. Any 
additional gluon emissions are not important for Xp > 0.8, since such processes lead to 
strange quarks with smaller Xp. 
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Figure |] and Table |T]f) show the cross sections for three Xp bins in the current region 
of the Breit frame for the full range, after removing the diffractive contribution with 
hmax < 2. The hatched bands represent uncertainties in the simulation of the 0-meson 
production by the MC models LEPTO, ARIADNE and HERWIG. The uncertainty due 
to As values between 0.2 and 0.3 is also included, such that the upper bounds of the 
hatched area for Xp < 0.8 correspond to LEPTO with As = 0.3, while the lower bounds 
of this area indicate the HERWIG predictions. For Xp > 0.8, the HERWIG prediction is 
between LEPTO with As = 0.2 and As = 0.3. The predicted cross sections of ARIADNE 
are always within the shaded bands. 

The MG uncertainties are small for Xp > 0.8. The predictions are shown with and without 
the contribution from the process of Fig. [^a). The measured cross section clearly requires 
a contribution from interactions with the strange sea. The MGs with the GTEQ5D or 
the MRST99(c-g) [Q (not shown) parton distribution functions reproduce the measured 
rate of 0 mesons. In these parameterisations, the strange sea is suppressed with respect 
to the non-strange sea, with the ratio ss/dd in the range 0.25 — 0.5, depending on x. The 
predictions correctly describe the results and thus conhrm the strange-quark suppression, 
even though the values of this data are signihcantly larger than the strange-quark 
mass. 


8 Conclusions 

Inclusive 0-meson cross sections have been measured in deep inelastic scattering for 10 < 
Q2 < 100 GeV2, 2 ■ 10“^ < x < lO'^, 1.7 < < 7 GeV and -1.7 < < 1.6. The 

MG predictions with a strangeness-suppression factor A* = 0.3 overestimate the measured 
cross sections. A smaller value of the strangeness-suppression factor, A^ = 0.22 ± 0.02, 
reduces the predicted cross sections and gives a good description of the total 0-meson cross 
section, as well as of the differential and x cross sections. However, Monte Garlo 

models based on Lund fragmentation fail to describe the and the Xp cross sections. 
The HERWIG simulation describes the measured cross section in the current region well, 
but predicts a smaller overall cross section than that measured; 0-meson production in 
the target region is underestimated by all MG models. 

The production of 0 mesons in the current region of the Breit frame has a signihcant 
contribution from the hard scattering of a virtual photon on the strange sea of the proton. 
The predictions for the rate of high-momentum 0 mesons with large values of the scaled 
momentum, Xp > 0.8, in the current region of the Breit frame have small uncertainties, 
since the 0 production in this region is dominated by 7 *s —>• s scattering. To reproduce 
the observed rate of 0 mesons at high Xp, the MG models require a signihcant contribution 
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from the strange sea of the proton. In this region, the measured cross section is correctly 
reproduced by these models when 7*5 —*• s scattering is included. These results constitute 
the hrst direct evidence for the existence of the strange sea in the proton ai x < 0.006. 
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a) 

range (GeV) 

da/dpif (nb/GeV) 

b) 

range 

da/dp'l' (nb) 


1.7-2.2 

0.433 ± 0.035l[J:[JJt 


-1.70- -0.80 

0.079 ± 0.009l){;|([(^ 


2.2-2.7 

0.237 ±0.019l°;°°2 


-0.80- -0.30 

0.175 ±0.0191°;°°^ 


2.7-3.2 

0.142 ± 0.016l|]:|^?^ 


-0.30-0.20 

0.199 ±0.0141°;°°° 


3.2 -4.0 

0.071 ± 0.009l|]:|j[jt 


0.20 - 0.70 

0.248 ± 0.022l°;°20 


4.0 - 7.0 

0.018 ± 0.0021°;°°^ 


0.70-1.60 

0.140 ±0.010l[););j(j 

c) 

range (GeV^) 

da/dQ"^ (nb/GeV^) 

d) 

range 

da/dxpftaiget) (nb) 


10-25 

0.02140 ± 0.001481°;™ 


0.0-0.5 

0.059 ± 0.009l°;°J° 


25-40 

0.00708 ± 0.00047l[J;S 


0.5- 1.0 

0.193 ±0.02ll[(;);^^ 


40-60 

0.00350 ± 0.00040l^;S 


1.0- 1.5 

0.112 ±0.0121°;°°° 


60-80 

0.00199 ± 0.00026l|];S 


1.5-2.0 

0.077 ± 0.010l)(;);[(^ 


80 - 100 

0.00079 ± 0.000211°;™ 


2.0-3.0 

0.049 ± 0.006l[););j(^ 

e) 

range 

da/dxpicniieYit) (nb) 

f) range(r7max > 2) 

dcr/cixp(current) (nb) 


0.00-0.30 

0.105 ± 0.022l°;°°? 


0.00-0.30 

0.105 ±0.0221°;°°® 


0.30 - 0.45 

0.263 ± 0.048l|];[J^^ 


0.30 - 0.45 

0.262 ± 0.048l[););^|? 


0.45 - 0.60 

0.203 ± 0.049l|];|J?^ 


0.45 - 0.60 

0.201 ± 0.048l)J;);^^ 


0.60-0.80 

0.135 ± 0.042l°;°26 


0.60-0.80 

0.121 ±0.037l[););| 


0.80-1.10 

0.090 ± 0.0141°;°}^ 


0.80-1.10 

0.072 ± 0.01ll[(;)5}^ 


Table 1: Differential <p-meson cross sections as functions of pif, rj'l’, and Xp. 
The statistical and asymmetric systematic uncertainties are shown separately. 
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10 < < 35 GeV^ 


range (all regions) 

da/dx (nb) 

0 .0002- 0.0006 

288.6 ± 33.01^5^8 

0.0006- 0.0010 

216.7 ±24.81^1 

0.0010- 0.0014 

110.0 ±8.5l|^ 

0.0014- 0.0018 

105.7 ± 17.31^1 

0.0018- 0.0030 

47.8 ± 5.3l^;? 

0.0030- 0.0050 

11.7 ±3.1+°;? 

Cnrrent region 


range 

da/dx (nb) 

0 .0002- 0.0006 

25.6 ± 23.41^;° 

0.0006- 0.0010 

57.4 ±9.9^^;? 

0.0010- 0.0014 

43.4 ± 9.4+°;° 

0.0014- 0.0018 

36.7 ±9.5+^;^ 

0.0018- 0.0030 

16.4±3.3l|)( 

0.0030- 0.0050 

4.5 ± 1.7+^;? 

Target region 


range 

da/dx (nb) 

0 .0002- 0.0006 

269.0 ± 33.4+?;® 

0.0006- 0.0010 

146.1 ± 16.61?;? 

0.0010- 0.0014 

78.5 ± 11.31?;? 

0.0014- 0.0018 

74.7±9.2l®;° 

0.0018- 0.0030 

32.4±6.6l?;? 

0.0030- 0.0050 

9.4 ±5.211® 


35 < < 100 GeV^ 


range (all regions) 

da/dx (nb) 

0.0008- 0.0015 

45.9 ± 13.011? 

0.0015- 0.0022 

39.0 ±4.61?;? 

0.0022- 0.0030 

36.0 ± 8.2l°;® 

0.0030- 0.0037 

22.0 ± 3.41?;? 

0.0037- 0.0060 

10.0 ± 1.81°;? 

Cnrrent region 


range 

da/dx (nb) 

0.0008- 0.0015 

13.0 ±5.111? 

0.0015- 0.0022 

16.4 ± 3.811? 

0.0022- 0.0030 

13.4 ± 3.811? 

0.0030- 0.0037 

11.5 ± 3.01?;? 

0.0037- 0.0060 

4.9± 1.2ll° 

Target region 


range 

da/dx (nb) 

0.0008- 0.0015 

34.3 ± 7.511? 

0.0015- 0.0022 

25.7 ±5.41?;? 

0.0022- 0.0030 

22.2 ± 6.2l?;? 

0.0030- 0.0037 

11.7±3.0l?;? 

0.0037- 0.0060 

7.2 ± 1.711? 


Table 2: Differential f-meson cross sections as a function ofx for two intervals in 
. The statistical and asymmetric systematic uncertainties are shown separately. 
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Figure 1: A schematic representation of different mechanisms for (j) production in 
inclusive DIS: a) a (p meson is produced from a strange quark after the interaction 
on the strange sea according to the QPM; b) a (p meson is produced from a strange 
quark emerging from the BGF process; c)-d) a (p meson is produced solely by the 
hadronisation process, independent of the flavour of the quark participating in the 
hard interaction. Additional sources for (p mesons are the hadronisation of strange 
quarks produced by higher-order gluon splittings, resonance decays, such as e) the 
Dg-meson decays; and f) diffractive cp-meson production. 


16 




































ZEUS 



> 

CD 


cn 

£i 

c 

CD 

> 

LU 



Figure 2: The invariant mass of the (p-meson candidates (points with statistical 
error bars) a) in the restricted kinematic regions pif > 1.7 GeV and —1.7 < rj^ < 
1.6; b) for the highest Xp value in the current region of the Breit frame, in addition 
to the cuts as in a). The solid lines show the results of the fit described in the text. 
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dc / d (nb / GeV^ ) do / d x (nb) do / dp^ (nb/ GeV) 
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• ZEUS 95-97 
LEPTO: 

I hadronisation 
s (s) from BGF 
I I s (s) from QPM/QCDC 
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Figure 3: Differential (j)-meson cross sections as functions of a) pif, h) rj^, c)-d) 
Xp and e) , compared to LEPTO, ARIADNE and HERWIG. The LEPTO and 
ARIADNE predictions are shown for Xg = 0.22. The data are also compared to 
contributions from LEPTO events with 0 mesons produced in hard interactions (s 
or s from BGP (light shaded area), from QPM/QGDG (dark shaded area)) and from 
events without strange guarks at the parton level (unshaded area). The full error 
bars include the systematic uncertainties, which are typically negligible compared to 
the statistical errors. 
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Figure 4: The inclusive sections as a function of x for two intervals, for 

the current, a)-b), and the target, c)-d), regions of the Breit frame compared to 
the HERWIG (dashed lines) and the LEPTO (solid lines) predictions with Xg = 
0.22. The LEPTO model shows separately the contributions from events with (f 
mesons produced in hard interactions (s or s from BGE (light shaded area), from 
QPM/QGDG (dark shaded area)) and from events without strange guarks at the 
parton level (unshaded area). 
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Figure 5: The cross sections for leading cj) mesons as a function of Xp in 

the current region of the Breit frame for r^max > 2. The hatched hands represent 
uncertainties in the simulation of the (p-meson production by Monte Carlo models, 
and include LEPTO (Xs = 0.2 —0.3^), ARIADNE and HERWIG. The upper bounds 
of the hatched area correspond to LEPTO with = 0.3, while the lower bounds 
of this area are defined by the LEPTO (Xg = 0.2) and HERWIG predictions (see 
text). 
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